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ABSTRACT 


The computation method developed in this thesis proceeds 
from the theory developed by Jones [Ref. 1]. fees eet drat 
eguations are first rewritten in matrix format. They are 
then organized into computaticnal Sequences that may be 
translated into computer programming language to calculate 


the aerodynamic parameters of wings of arbitrary planforn. 
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I. BACKGROUND 


The cverall research of which this thesis forms a part 
has a three-fold purpose. Petast ly, it is intended to 
instruct the student about the fundemental equations that 
govern the aerodynamics of wings in incompressible flow. 
Secondly, it teaches the student how to convert the funde- 
mental equations into a form Suitable for numerical calcula- 
tion on the digital computer and to carry out the complex 
programming required for this  furpose. Thirdly, once the 
first two objectives have been accomplished, the computer 
program itself provides a useful pedaygoyical and design tool 
for illustrating the effects of various wing design parame- 
ters on the final aerodynamic performance of a wing. 

Lcedr. J. Le. Parks made a creditable start toward the 
first objective in his thesis [Ref. 2]. However, in the 
limited time available to him, he was unable to make signif- 
icant progress toward the second or third objectives. 

It was evident early that Parks' computational technigue 
had some flaws. The expectation was that these would be 
readily discovered and soon ccrrected and that objectives 
two and three would be rather guickly and easily completed. 
It has not worked out that way! 

As the investigation actually developed, a whole series 


of obstacles were encountered and each was in turn eventu- 


ally overcome. But the process was neither quick nor €asy. 
In every case, 1t was found that the seeming obstacle was 
really attributable to some ccnceptual error. Once the 


error was found, the obstacle disappeared and the investi- 
gator gained a deeper insight. 
Among the obstacles that were encountered and eventually 


overcome were the following; 


1) Confusion about the proper type ol weenpies Grecia 
grid to use, whether staggered or unstaggered, how 
fine, and so on. We have now settled, for good 
reasons, on an8 by 8 unSstaggered grid over the 
semi-Span. 

2) Confusion about how to deal numerically with the 
singularity that occurs in the governing integral 
equations. One option is to try to avoid the 
problem by staggering the grid of the field points 
with respect to the grid of the control pcints. 
Instead of this, we have resolved the indeterminacy 
by rigorous analysis ard as one consequerce, are 
able to employ a simple unstaggered grid. 

3) Confusion about how to evaluate the partial 
derivatives of the circulation function, whether 
analytically or by finite difference formulas. it 
was found that analytical differentiation iS incor- 
rect and that finite differences must be used. 

4) Confusion about the validity of representing the 
Circulation by a Fourieé€r series. While this proce- 
dure is widely advocated in the technical litera- 
ture, we found that a much simpler and clearer 
formulation can be obtained otherwise. 

5) Confusion about the boundary conditions. Special 
conditions apply to the leading and trailing edges, 
the wing tip and at midspan. These quite complex 
boundary conditions have now been fully and rigor- 
ously analyzed and incorporated consistently into 
oul farmulaticn of the froblen. 

Jt is evident from the foregoing discussion that this 
res23rch has amounted to a majcr education in baSic aerody- 
Namie and it numerical methods. In these respects it has 
reen a deeply rewarding experience. However, in view of the 


for: joing obstacles and problems, the time schedule has of 


course been greatly delayed frem that which was initially 
anticipated. Thus there are nce final numerical results at 
this particular stage of the investigation and there is no 
more time available to the present reSearch team. Hence, 
this final aspect will have to be completed by some subkseg- 
yent investigators. Nevertheless, what the present effort 
haS produced to begueath to any subseyguent worker is a 
sophisticated and refined numerical method. This method has 
now evolved to the point that it can be confidently expected 


to produce a reasonable and accurate final result. 


II. INTROLUCTION 

This thesis is the development of a computational method 
tor determining the aerodynamic parmeters of wings ihn incon- 
Eressible flow, from the theory developed by Jones [Ref. 1]. 
In conjunction with Jones" theory this analysis has solved 
the problems reguiring further attention pointed out by 
Parks [Ref. 2]. Primarily a new approach, defining the 
circulation directly instead of uSing a Fourier series 
representation, is employed. This along with a direae 
method of dealing with the singularities that are encoun- 
tered in the joverning integral eyuations when a_ control 
point and field point coincide, allows the computational 
grils for those points to be identical. This analysicme. 
intenled to complement that of Jones by taking his final 
theoretical e,uations and, where applicable, presenting then 
Ln ier I ees tee The eguations are then oryanized into a 
comrutational sequence which can be further translated into 
programming language. It 1s expected that a program to 
produce those wing parameters described by Jones will fcllow 


from this development. 
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Tif. SATRIX NOTATION 

Though equations using subscripted variables are a 
brinary method of portraying mathematical relations, it may 
92 easier to visualize what is happeniny in those equations 
waem they cre written in matrix format. In this chavter 
those eaqvations of reference 1 which Lend themselve to this 
type of notation are rewritten in matrix format. The 
Subscripts (indices) i,j,k are described for the wing in the 
Ge plane <zs shown in figure 3.1. 

Tt will be convenient at many points throughout this 
presentation to transform two dimensional matrices into 
vectors and conversly vectors into two dimensional matrices. 


For doing this the following sytolic relations are defined; 


>| | ie (3.1) 
@xg  ? oxV (4x! 


K 


Figizre 3.7 as it is drawn represents atwo dimenSional 
Matrix with idices (J,I) and the index K Shows the relation 
of a column vector that has been transformed into a matrix. 


The inverse operation can also te deduced from this figure. 


A. CIRCULATION AND PRESSURE REIATIONS 


The following equations are the general circulation 
derivatives and pressure relaticns. Those eguations taken 
From reference 1 are noted by their ejuation number preceded 


by tuc letter J. Therefore (J6.8) becomes; 


yO ; l 

a \ > (Pf 

ia. . g p ied) 
J 5x8 64x! 

Slmilaniy for (J6.9) 
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OE Rene ates» ee i ee we es St Cee eee ee er eip etme i Oe See Se a | 





Figure 3.1 Index notation. 


(35), le}? ae} a 


64x | 
Along the leading and trailing edges equations (J6. 10) 


through (J6.13) respectively are 


as 
oes 4 ~ 1B} (3.4) 





: Gx) 
OPN. 
(S6),47 Spr} _ 
, FX! 
lenin 
(S92), 5° Tow3 a 
26 - 
(S5),5 2 oer} io 
5x 
The ,-ressure distribution (J6.14) is 
ACp)4.i =P(4,i) = | P (3.8) 
GX 


At the leading and trailing edges (J6.15) and (J6.16) are 


ACp)j= AG)> RS (3.9) 
ENC») 74 ae LAS (3. 10) 


Generally the pressure eguation (J6.22) can be written; 


[P= al elt |Pame am 


or solving for the circulation derivative in the chordwise 
(J) direction 


J) 


]=SPlLP |p 


8x Gx6 gxy 8x (3.12) 





Be. CALCULATION OF DERIVATIVES 


1. Chordwise Derivatives 


Eyuations (J7.9), (J7.1€) and (J7.14) are written in 


the index notation as follows 


5G, N=a LPG $0 2) (3. 13) 
Rada lalGelteC Gi) i237 Bey 
(4.8) <a 3G, DB, 82] (3.15) 


This helps to visualize the eyuation in the matrix notation 


below 


[fe] =D. Jr oe 


Sx& 5x 


ls 


Note that this equation involves the transpose of both theft 
matrix and the fg matrix. A description of the matrix Da may 


be, found in figures... 


an ee me SR GET Lee ES SOS) Ge EE NS EP eS 


1 $lelelefololo! 
zfeole [olejole fe 
+/ 
3 o|sle 


| 

, 

| 

| t 

_ 1 [ele 

lo. J: a5 SE Zo 

! ofz le 
| 

















EEE 
~ t! 

O 

eo fojololololaly. 
Figure 3.2 Matrix Da. 


Equations 3.13 through 3.15 can also be written ina 


Single index notation. 


} 
(3 (k)= ag lP Ck) +3 r(he8) | 1Sk$8 (3.17) 
lo (k) = aa lal (k-8)* 20(k+ 8) | 9<k<56 (3. 18) 
lo (h)= aol 3 (h-8)+ 2P(k) J S7SkS64 (3.19) 
In order to get adjacent points in the chordwise direction, 
the integer 8 must ke added to or subtracted from the k 


index describing the foint as can be seen in figure 3.1 The 


above equations now can be written in matrix notation 


bes = [D2 Ji} a 


64x| 64x64 GYx! 
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where the general configuration and description of the 
elenents in the matrix Da~ ace shown in figure 3.7 and 
table 1 respectively. Notice from the figure that all of 
the nonzero elements of the matrix Da™ lie along just three 
Perncipal diayonals. 

Since in the basic lift case the pressure at the 
leading edge is finite, an additional term appears in the 
equation for the circulation derivatives. That term 
involves the second derivative cf [. Since lg is zero, equa- 


tion (J6.28) for the pressure at the leading edge is, 


_ si 

LAS =a LE) {les = 
&x | Sx | 

or solving ae: 


UToes5 = #[¢ J UP (3. 22) 


& x8 Gx! 
By nc the elements of both {fjz,} and {P,} and placing 
them in the first cclumn of the partitioned matrices with 


the same name thusly, 


[eee ] | {fo0.$![ 0 a eae 


Sx y Sx! 


| 
={{P5'1 0 | 
[RJ [{PJilo]| a 
5X§ sx! Rx] 
the vector eyuation can then be written in matrix form 


Tou J= % a) ‘ [PJ ih 


Now the [yg eyuations, (J7.10), (37.16) and (J7.14) are; 
L\E 
7 leer Bi )* Toy! E 206 PCy A) | (3. 26) 


15 


[s(y,4) Eon Pai My aeal Geel] 1=2, 3g 7! 3.01290) 
Io (4,8)=a0l3l(y,7)* 4 r(4,8)]J (3.28) 


for which the matrix notation is; 
D r T T 
F | loek| ni) = | D Ie] 
7 z 2 (3.29) 
ox sxx &xKS &xXS 
or solving for circulation, 


[r] = [D.y (42h) * fhe] ) (3. 30) 


& XS 


with the contents of Db shown in figure 3.3 





Figure 3.3 Matrix Db. 


Notice that the element of Db are identical to those of Da 
except cor the first two elements of the first row, which 


are 6 and 0 vice 1 and 1/3. 
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2. Spanwise Derivatives 


For the derivatives in the spanwise(@) direction a 
set of relations similar to those for fy can be shown. Froa 


(J7.17), (37.18) and (J7.23) respectively 
o(1,4)2alFC,4)+39(2,4)] (3.31) 
ee ee ee yee a 
[o(4.4)=20[% MG Wi )e2P (get, 4) ] pees see) (2 32) 
BS, 4)=gelz1 (6, «) -20(7,4)+20(3,4) | (aeyce) 


So that in matrix notation the equation is, 


fe J= CE JIN ae 


and the E matrix is defined in figure 3.4. 





Figure 3.4 Matrix E. 


BI 


In a manner similar to the [% equations, eguations 


3.31 through 3.33 can be writter in single Side woe nod 


[o(k)=40l P(h) + 3 CCk+1) | K=1,9,17,002,57 (3-35) 
[e(k)=z5] 31 (k-1)* 2 (hel) ) k=2--7,10--15,...,58--63 (3. 36) 
[ol k)=g@ [21 (k-2)-20(k-1)* 3 7(k)) k=8,16,24,...,68 (Secu 


These eguations when viewed in the vector form are much 


neater. 
h(E HS (3-38) 
eee “CHKE CE O¢xI 


The yeneral configuration of the E* matrix along with the 
description of its elements is shown in figure 3.8 and table 
Ze Notice from the figure that all of the nonzero elements 
Ob) Welle lee lie along just four principal diayonals. 
Storage can be minimized by storing [Da*] and [Db*]} as 64 by 
3 matrices and [E*] as a 64 by 4 matrix, capitalizing on the 
large numker of zeros in each. 

It is useful to consider the hypothetical case in 
which the circulation function is assigned unit value at an 
arbitrary foint Q (the ‘circulation point") with matrix 
indices JQ and IQ or equivalent vector index KQ, and is 
assityned zero values at all other points of the calculation 
Cigars 

The above hypothetical distribution of [| implies 
that of the many field points (denoted by index kK) there 
will ke certain ones adjacent te the circulation point Q at 
which tke derivatives [y and [.take on nonzero values. in 
the simplest and most typical case where Q does not lie 
adjacent to any of the boundaries, there will be two field 


points, at locations K=(KQ-8) and K=(KQ+8), at which [yg takes 
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on the values +1/2A@ and -1/2O¢. fhicew toe me here: will be 
two other points, at locations K=(KQ-1) and K=(KQ+t1), at 
which f, takes on the values +1/246 and -1/2Ae. This situ- 


ation iS illustrated by the schematic diagram in figure 3.5. 








<<: SSS OO eS ee eee 


we (ke-): 330 







¢/ 
lo \k@- &) = 24g k= KQes 


(Kaye 2, kas)* 350 





K~KQ¢ti 
ae a) ime 











Figure 3.5 Circulation point. 


More complicated situations arise if the point Q 
lies immediately adjacent to a boundary. In this case the 
humber of field points having ncnzero values of [gand [6 may 
range from three to five. The details can be seen fron 
Careful study of the matrices Da’, Db™ and &*. 

Along the trailing edge a set of relations for the 
Gieculation function and its Spanwise derivative is 
reguired. From a Similar analysis to that used for the 
Spanvise derivative over the wing, the following set of 


equation for the trailing edge is obtained. 
q | 
Jor (= gel) sR (2)] (3. 39) 


fe =a EGG-O teh 5-0] j=2,3pee047 (3440) 


V3 


as )=hel = Pie ee ce (3.41) 


This set of eqguations can be written aS a single matrix 


equation 
{Ter ; LE M3 (Sa 


This E matrix is the same as the one described™ ecarlien a 


figure 3.4. From equation (J7.15) at the trailing edye 


ot oa) = eee , oF 
p= ae g.7)+ 078) (3-43) 
Rewritten as a vector equation 
C(49) r(57) 
(03+ alriol Bl 
ihe ~ 26( (56) 26 ( F(¢4) (3.44) 
5x | Sx | 5x! 


Notice that of the 64 elements in the [” vector only the last 


16 are used. 


Ce INTEGRAL EQUATION FOR WING SLOPE 
|. JEL Stee any Slepem@ tune CLeh 


The integral eyuation (J7.24), 1s restated here as 


ejguation 3.45 


iT 
a 8, +6 IE ~ fe), (nee) S 
W/o = ITA Ee Pe ine |, ieee ok sepa 5 sine I, ddde (3-45) 
OSS “ 


By defining the following two relations 


_ apa f f 2] 
T( ke, k) = een | i 53 |SING (3.46) 


and 
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S (ke, k) = poe | CCP. Tee) |S sing (3.47) 


the integral equation can be written in index notation as 


C4 «64 
\op (kp) = (> i: (7 kp,k) D*( ky ka) - SC kp, RIE 7) CO) (3.43) 
Roz/ keel 


The matrix D“ represents either Da* or Db™ whichever is 
appropriate. The elements of Da™ and Db* are shown in takle 
1 and those of E”™ in table 2. 

To make the distinction between additional and tasic 


lift the ae twc matrices are defined 


AN Ue, he) =, Y (Ty A) Oe (Rka)-S(ha KIER ke)) (3.49) 


B (kha) = 2 (Cte k)D%(K ke)-S(hek)ERAa)) (3.50 


The above caculation for [A']} and [B'] is simpler than the 
64 by 64 multiplications indicated. If one uses tables 1 


and 2 the products are simplified to, 


6¢ 


T (ke, 2) k) Dg CR, ka) = T (kp, ka-s) 0 (ke-S, ka) 
. + T( kp, hq) DE Cke, ke) 
+ T (kp ka+8) Di (kas, ka) (3.51) 
C4 
DT (hy, k) E(k, ke) = Thy, ho) DE (has, a) 


+ T(kp, hq) Dé (ka, ko) 

tT kp, Ra&) OX (korh ko) (3-82) 
GH 
Y. SCkp RE“ (k ke) = S Che, ha-t) E*(ke-t, ke) 
hel 

t S (hp, ha) E*(ke,ke) 
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+ S(ke ho+t)E(Ro+l ko) 
¢ Se ko+2) E (hax, ko} (3. 53) 


The calculation of [A‘*] and [B*? then is 


A‘(ke ho) = T( ke, ha-8) Di(ka-s, ke) 

+ T (kp kq)D2 (ka, ko) 

+T (RP ka+s) Oe (Rars, ke) 

aS (kp kg-l) E’(ko-t, ke) 

-S (ke ka)E* (ka, ko) 

-§ (kp, katl ) E"(ka+l,ka) 

-S CRP ka@e2) E*(Rat2, ko) (3. 54) 
Py'( ke ko) = TC ke ba-&) Deka -s, ke) 

+T (ke ke) O¢ (kha, ko) 

+T (Re kats) Ob kOe ko) 

-S (ke ke-l) E’(ka-t, ka) 

-S (ke, ka) Eke, ko) 

-SCkp, kat!) E*CkG+t ka) 

-S(he ko-2) E*(h@+2, ke) (3.55) 


The various factors that occur in the above equations take 


values if and only “ae 
and KQ are first assigned 


on definite numerical the indices KP 


definite numerical values. when 


equations define the value of 
Single element of either [A'] cr [B"}. 


this is done the above a 


The matrices [A'] 


and [{B'j] can be stored as two 64 Ly 64 arrays. So the 


matrix equations for the first wing slope function are, for 


the additional and basic lift cases, respectively, 


we} = [A J irs (3-56) 


64x 61x64 64 xI 


ae 


64x64 64x 


be = [BJU bs 


2. Second Wing Slope Function 


The integral equation for the second wing’ slope 


function (J7.25) is 


% 
_ 7) | G, 4 G2 | [, 
7 = ae | (y-%) (qenjJ ler de (3. 58) 


This equation requires only the derivative of the circula- 
tion function in the © direction along the trailing edge. 
Therefore a relation involving the same E matrix as before 


can be utilized. By defining a function 
(3259) 


HC ke, 4) = - tk re Kel + Gs) 


equation 3.58 can be written in index notation as 
§ 
wehbe) =) He) Der (y) (3.60) 
FF 


In Matrix notation it becomes 


S23 = LH] $ 3 (3. 61) 


64x! 6YxS = §xI 
Equation 3.42 can be substituted for fer and the equation now 
becomes 
ff 
we} =[4 JE] £03 (25 62) 
64x | 


xf xi 


Then substituting equation 3.44 for [+ gives 


ws jn 
(6) = [HIE Ct rad “Ble bs 


67x! LH Xxs 


For convenience an inns matrix is es ai as f£cllows, 


Zo 


[a]: [HCE 2.0 


CG4Yx¥ §xZ 
Now another matrix can be defined which shall be termed the 
Vomatrix. Matrix Vis a partitioned matrix, which is 


divided into three parts as shown in figure 3.6. 





69x45 


Figure 3.6 Matrix Vv. 


Therefore the final form of the integral eguation is, 


hwo §=[v] IF 2.65 


64x! 64x64 67x! 
3% 


rr 


inal Wing) slopesrunerien 


The two relations for Wr' and Wp" can now be super- 
imposed to give the final equation for the wing slope func- 
ELOne es FO mel? 42) 


os Th 
vip § iis ; \/o * Swe (3. 66) 
Cr ox : 64x] 


After substituting the relations for Wp' and Wp" and simpi- 


fying, the eyuations become 
ae | 
i wos = (LA J +[v]) £73 (306m 
64x! 64x64 64 xhy4 6YX] 
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We 3 > 3+ ([B]- vq} 3 


(3.68) 
CYx | 64x64 fx bY GY x 


By defining two new matrices [A_ 


- and [BJ] as follows, 
Al= [a‘}-[v] 


64x64 


64x64 64x64 (3. 69) 
| B fe B’ + 
V (3.70) 
64x64 64x64 64x64 
the final 


wing Slope functions can be written as 
(we = [A] (FJ a7 
bYx | 64x64 6x I 


for the additional lift, or as 


{we} =(B] {3 


(se 2} 
G4¥x| 64x64 by xl 


HoE the basic lift. The additional lift eyguation 


can be 
inverted as follows to be used in later calculations 
; 7 -| 
| | A | Wp § (3.73) 
Ca] (xo CYx] 
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IV. SUMMARY OF KEY RELATIONSHIPS 


Ae PRELININARY CALCUIATIONS 


The initial set of calculations to be performed by the 
computer are for the values of those constants which need 
only be calculated once. The following equations use only 


the input parameters AR,2A,7 and the indices i and j. 





Oe Vs (4.1) 

We = 6 (4.2) 
,_ dAdAde 

C = 217m (4.3) 
Ae 

e- 2TTAR (4.4) 


The following equations involve vectors whose indices i or j 


go from 1 to 8. 


D(x) = AGlA-$) (4.5) 
ea) - Ae: > (4.6) 
ames ) = #(1-cos(gr))} (4.7) 
1 (7)= Cos(ey)) (4.8) 
Clg )= mR (i+ r)-22-yG)) (4.9) 
ee) = in IEE (i) (4. 19) 


Sal 


Diz)= O (4.11) 


Dr (7) = 1T (4.12) 
O,(i)= O (4. 13) 
nae wee (4. 14) 


Indices i and j can be obtained from the index k through the 


two following integer equations 


i 


wl (R/g) +1 (integer division) (4. 15) 


i 


4 k -3 (4-1) (4.165 


Indices ip and jp can be obtained from kp similarly. 


B. CALCULATION OF THE A AND B SATRICES 


1. The At and Bt 


atrices 


Since the matrices [A} and [B] are the sum of the 
[V} and ;A‘'j or [B']J matrices, each of the matrices will be 
computed separately and then brought together to form [A] 
and [BJ]. Recall the equations 


PA'J= [7 Jfor}-[s}l "J (4.17) 


8 J=[T JloJ-[s]le" 2.10) 


Where [TJ], [5], [D*] and [E*] have been previously defined. 
It 1s convenient to make the following variable 


definitions 


SZ 


A”) = (1- Te) (4. 19) 
Ap =a (1+ He) (4.20) 
NK = \(An)+C5- p&p (ie) 


= VAX" + An" (4. 22) 
/ 2 
i a\x “AN, (ee 


F,= AX ~BN, 4 (4. 24) 

E = AX- SI A (4. 25) 
T (kek) =e ‘(as » Be) sin OY) (4. 26) 

SCke, k) = C (45 y 4%) SP! ci Bea) (4. 27) 

(4) 

T (Rr, ke)= C \ 43 / SING (yp) (4. 28) 

Anr\C ly) 

S (kp, | 2p) aC (3) SF sin OU) (4. 29) 


A check sheull be made to see whether the control point and 
Peta point coincide. If they do, then equations 4.28 and 
4.29 shail be used instead of 4.26 and 4.27 The matrices 
[A']} and [}3'] from eguations 3.54 and 3.55 are 


A (ke ko) = T (ke, ka-8) Oo(ka-& ho) 
+ T (ke, ko) D&C ko, ko) 
+T (Kp kat) Del kes, ka) 
-S (kp ko-t) E*(ka-l, ko} 
-S( ke R@)E* ko, ka) 
-$(Rp, kot!) E*(ka tl, ka) 
~S( ke ka+2) E*(ka+2, ka) et 
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B'( kp ka) = T( ke, kq-&) De(ka-8, ha) 
+T (ke ko) DECRG, ka) 
+T (ke ka+3) Di (Rots, ka) 
-S (kp, ka-!) E*(k@-l ke) 
-S( kp RaQIE*( ka, ka) 
-S( kp kot!) E*(kat], ka) 
SC ke RQ+2) E*(kat2, ko ) (4. 31) 


The various factors that occur in the above equations take 
on definite numerical values if and only if the indices KP 
and KQ are first assigned definite numerical values. When 
this is done the above equations define the value of a 
Sinyle elenent of either [A‘'] cr [B']. The matrices [A* } 


and ; B'j can be stored as two 64 by 64 arrays. 


Recall figure 3.6 with the associated eguation 3.64 


repeated here to provide clarity. 


eas [f° o 1% [a :"#( a] (4. 32) 


64x64 
LQ] = rut e} (4. 33) 
6Yx¥ 64x5 sxe 


In addition to selected definitions above these additional 


definitions are made 








DK = NA) eee (4. 34) 
is Vaxe+an? (4. 35) 

= VAXP+An, (4. 36) 
G,-1- os (4. 37) 
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= i L\ Xr 
cr Fo (45373) 


H( ke ¢)= 0’ an* ahs) (4.39) 
4 (ke 4) =C (CH } (4.40) 


A cieck Sheullt be made to see whether the spanwise coordi- 
Mat?: coincice. If they do then eyguation 4.49 shall be used 
instead of 4.39 Then [QJ] can be computed 

x 


QC ke, 4) =). H (hp, 4) E(g,4) (4.41) 


Now the [V]} matrix can be constructed as was shown in figure 


0 « Since the first 48 columns of V are zeros, only the 
nonzero elements of V need be stored. These amount toa 64 


by 16 matrix. 
3. Zhe A and B Matrices 


Finally these matrices are superimposed to give the 


Aani3B matrices. 


A (hp, ko) = A (kp, ko) + V¢ ke, he) (4. 42) 
B(kp ka) =B (kp ho) +Vhp, ka) (4. 43) 


This gives the complete A and B matrices, which will be 


utilized in the following calculations. 


So ADDITIONAL LIFT 


Pommtne specific additional lift, as explained in more 
detail by Jones [Ref. 1], setting the normalized slope to 


-1, gives 


32 


313 ‘ | 2 rE (4. 44) 


64x G¥x6Y 64x} 


This circulation can then be transformed by the operation 


described earlier to get 


(3 9 £0] _ 


64x! 


After being transposed it can re used in equation 3.16 to 


get the chordwise derivative 


a) = (oJ iry oi 


SKS 
The corresponding derivative in the © direction is not 
generally needed and is not shown here, but has been stated 
earlier as eyuation 3.34. Now es transposed and subksti- 


tuted into equation 3.11 to get the pressure distribution. 
ic Geel XN 4 
je J=2 eI (] End | (4.47) 
XS xy SxS SXS 
From this pressure distribution one can calculate the aero- 
lynanmic certer and the slope of the wing lift curve. 
ar Pxoleet ir. 


In the case of basic lift, the pressure distribution is 


initially —~rescribed. From thie dvStribution onewevtains. 


TeJ= 4 3 0s c JL p ]['sms (4. 48) 


Fx §FxsZ 


and additionally along the leading edge 


| ope Bf c|] 2 | (4.49) 


ake: SXF Ex 


S10 


where only the first column of each of the matrices [fy,,] and 
[P,] is nonzero. Then [fog] and [fj] are transposed and 


substituted into equation 3.30 to jet 
- at 1 /ad a ci 
|| ] : [De | 4 [ewe J*L fe J eee 
SS aed SKS SxS 

The transpose of [f']’ is the basic circulation which after 


being transformed into a 64 by 1 vector is utilized in the 


following eguation to get the wing slope function. 


Swe3= [6] {P3 (4-51) 


61x] 64x64 64«I 


The wing slope functicn can then be used to find the varia- 
tion in the height of the mean wing surface above the refer- 


ence xy plane. 


ay 


V. CONCLUSIONS 

From the equations developed and listed in the summary 
of key relationships of the previous chapter it can be seen 
that they may be translated into a computer program which 
calculates the wing slope function and the lift distribution 
of any designated wing. A further straightforward extension 
of tie analysis (not given in this thesis) then fixes the 
various aelodynamic parameters of interest Such as the slope 
of the wing lift curve, the location of the aerodynanic 
center, the pitching moment ccefficient about the aerody- 
Namic center, the induced drag and the variation in the 
height of the mean wing surface above or below the reference 
XV plane. 
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